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1. INTRODUCTION 


This document is the final report of a research project 
concerned with the optimal design of helicopter rotor blades. The 
report contains three main parts: 

1. a discussion of the reasons for which the research 
was undertaken; 

2. a summary of project accomplishments, presented in 
the form of a list of optimization problems which 
have been solved and a list and brief description 
of findings related to optimization of rotor 
blades ; 

3. the doctoral thesis of Timothy Ko, which contains 
many details of the computations performed during 
the project. 

2. IMPORTANCE OF THE RESEARCH 

The design of helicopter rotor blades involves not only 
considerations of strength, survivability, fatigue, and cost, but 
also requires that blade natural frequencies be significantly 
separated from the fundamental aerodynamic forcing frequencies 
(e.g. Ref. 1) . A proper placement of blade frequencies is a 
difficult task for several reasons. First, there are many forc- 
ing frequencies (at all integer-multiples of the rotor RPM) which 
occur at rather closely-spaced intervals. For example, 5/rev and 
6/rev are less than 20% apart. Second, the rotor RPM may vary 
over a significant range throughout the flight envelope, thus re- 
ducing even further the area of acceptable natural frequencies. 
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Third, the natural modes of the rotor blade are often coupled 
because of pitch angle, blade twist, offset between the mass 
center and elastic axis, and large aerodynamic damping. These 
couplings complicate the calculation of natural frequencies. In 
fact, the dependence on pitch angle makes frequencies a function 
of loading condition, since loading affects collective pitch. 
Fourth, the centrifugal stiffness often dominates the lower 
modes, making it difficult to alter frequencies by simple changes 
in stiffness. 

In the early stages of the development of the helicopter, it 
was believed that helicopter vibrations could be reduced (and 
even eliminated) by the correct choice of structural coupling and 
mass stiffness distributions. However, it is easy to imagine how 
difficult it is to find just the proper parameters such that the 
desired natural frequencies can be obtained. The difficulties in 
placement of natural frequencies have led, in many cases', to 
preliminary designs which ignore frequency placement. Then, 
after the structure is "finalized" (either on paper or in a 
prototype blade), the frequencies are calculated (or measured) 
and final adjustments made. Reference (2) describes the develop- 
ment of the XH-17 helicopter in which a 300-lb weight was added 
to each blade in order to change the spanwise and chordwise mass 
distir ibution and thereby move the first flapwise frequency away 
fromn 3/rev. However, these types of alterations are detrimental 
to blade wight, aircraft development time, and blade cost. In 
addition, corrections usually are not satisfactory, and the heli- 
copter is often left with a noticeable vibration problem. 
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The state-of-the-art in helicopter technology is now to the 
point, however, that it should be possible to correctly place 
rotor frequencies during preliminary design stages. There are 
several reasons for this. First, helicopter rotor blades for 
both main rotors and tai: rotors are now being fabricated from 
composite materials (Refs. 3 and 4). This implies that the 
designer can choose, with limited restrictions, the exact El 
distribution desired. Furthermore, the lightness of composite 
blades for the main rotor usually necessitates the addition of 
weight to give sufficient autorotationa 1 blade inertia. Thus, 
there is a considerable amount of flexibility as to how this 
weight may be distributed. Second, the methods of structural 
optimization and parameter identification are now refined to the 
point where they can be efficiently applied to the blade struc- 
ture. Some elementary techniques have already been used for the 
design of rotor fuselages (Ref. 5). It follows that the time is 
right for the use of structural optimization in helicopter blade 
design. Some work on this is already under development, and, 
although not published, some companies are already experimenting 
with the optimum way to add weight ot an existing blade in order 
to improve vibrations. 

The purpose of the research project described in this report 
was to investigate the possibilities (as well as the limitations) 
of tailoring blade mass and stiffness distributions to give an 
optimum blade design in terms of weight, inertia, and dynamic 
characteristics. The work has focused on configurations that are 



simple enough to yield clear, fundamental insights into the 
structural mechanism but which are sufficiently complex to result 
in a realistic result for an optimum rotor blade. 

3.0 SUMMARY OF THE PROJECT 

3.1 OPTIMIZATION PROBLEMS WHICH WERE SOLVED 

The basic structure optimized was a beam free at one end and 
supported at the other. Various support conditions and con- 
straints on natural freqencies were used. The behavior of the 
beam was computed by using a 10-element finite-element model. 
Quantities associated with the finite-element model, such as the 
thickness or area moment of inertia of each element, served as 
design variables in the the optimization procedure. A typical 
formulation of an optimization problem was 

Find the flange and wall thicknesses of a box-beam 
cross-section (three variables per finite element) 
which minimize the weight of the beam, while main- 
taining the first natural frequency within a "win- 
dow" (e.g., 2.4 < p^ < 3.0 per rev). 

All optimization problems were solved with the CONMIN com- 
puter program [6]. CONMIN is based on the mathematical nonlinear- 
programming method of feasible directions. 

The list of problems solved follows. 

Case 1 . Cantilever beam 
Rotating: no 

Objective function: weight 
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i Design variables: area moments of inertia 

r 

Boundary condition ( s ) at root: fixed 

| Frequency Constraints: first flapping specified through equality 

* 

j constraint 

Autorotation constraint: no 
* Stress constraint: no 

References: First Semi-Annual Report pp. 17-18, Thesis, pp. 18-19 


Case 2. Cantilever beam with tip mass 
Rotating: no 

Objective function: weight 

Design variables: cross-sectional areas 

Boundary condition (s) at root: fixed 

Frequency Constraints; lower bound on first flapping 

Autorotation constraint: no 

Stress constraint: no 

References: First Semi-Annual Report pp. 17-19, Thesis, pp. 18-21 


Case 3 . Wind-turbine blade 
Rotating: yes 

Objective function: weight 

Design variables: area moments of inertia, lumped weights 
Boundary condition(s) at root: fixed 

Frequency Constraints: windows on first and second flapping 

Autorotation constraint: yes 
Stress constraint: no 

References: First Semi-Annual Report pp. 21-24, Thesis, pp. 34-37 
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Case 4. Hingeless rotor-blade 
Rotating: yes 

Objective function: weight 

Design variables: area moments of inertia, lumped weights 
Boundary condition (s) at root: fixed 

Frequency Constraints: windows on first and second flapping 

Autorotation constraint: yes 
Stress constraint: no 

References: First Semi-Annual Report pp. 23-27, Thesis, pp. 38-41 

Case 5. Cantilever beam with two frequency constraints 
Rotating: no 

Objective function: weight 

Design variables: area moments of inertia 

Boundary condition (s) at <:oot: fixed 

Frequency Constraints: windows on first and second flapping 

Autorotation constraint: no 
Stress constraint: no 

References: Second Semi-Annual Report pp.7-16, Thesis, pp. 23-33 

Case 6. Cantilever beam (similar to Case 5, except for three 
rather than two frequency constraints) 

Rotating: no 

Objective function: weight 

Design variables: area moments of inertia 

Boandary condition (s) at root: fixed 

Frequency Constraints: windows on first, second and third flapping 

Autorotation constraint: no 
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Stress constraint: no 


References: Second Semi-Annual Report pp. 17-18 

Case 7. Cantilever beam (similar to Case 6, except for addition 
of lumped weights as design variables.) 

Rotating: no 

Objective function: weight 

Design variables: area moments of inertia, lumped weights 
Boundary condition (s) at root: fixed 

Frequency Constraints: windows on first, second and third flapping 

Autorotation constraint: no 
Stress constraint: no 

References: Second Semi-Annual Report pp. 17-21 

Case 8. Cantilever beam (similar to Case 7, except for addition 
of autorotation constraint) 

Rotating: no 

Objective function: weight 

Design variables: area moments of inertia, lumped weights 
Boundary condition(s) at root: fixed 

Frequency Constraints: windows on first, second and third flapping 

Autorotation constraint: yes (constraint applied to mass moment 
of inertia of whole beam) 

Stress constraint: no 

References: Second Semi-Annual Report pp. 17-22 

Case ^Cantilever beam (similar to Case 5, except beam is 
rotating) 

Rotating: yes 
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Objective function: weight 

Design variables: area moments of inertia, lumped weights 
Boundary condition (s) at root: fixed 

Frequency Constraints: windows on first, second and third flapping 

Autorotation constraint: yes (constraint applied to mass moment 
of inertia of whole beam) 

Stress constraint: no 

References: Second Semi-Annual Report pp. 17-24 

Case 10. Teetering rotor 
Rotating: yes 

Objective function: initially the weighted sum of squares of 
differences in frequencies; after a feasible design is 
found, the objective is changed to the weight. 

Design variables: flange thicknesses, lumped weights 
Boundary condition (s) at root: fixed 

Frequency Constraints: windows on first, second and third collective 

flapping 

Autorotation constraint: yes 
Stress constraint: no 

References: Third Semi-Annual Report pp.5-6. Thesis, pp. 42-49 

Case 11. Teetering rotor (similar to Case 10, except cyclic 
flapping modes considered, instead of collective modes) 

Rotating: yes 

Objective function: initially the weighted sum of squares of 
differences in frequencies; after a feasible design is 
found, the objective is changed to the weight. 

9 

Design variables: flange thicknesses, lumped weights 
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Boundary condition(s) at root: pinned 

Frequency Constraints: windows on first, second and third cyclic 

flapping 

Autorotation constraint: yes 

Stress constraint: no 

References: Third Semi-Annual Report pp.7. Thesis, pp. 50-51 

Case 12. Teetering rotor (similar to Cases 10 and 11, except that 
both cyclic and collective flapping modes are considered) 

Rotating: yes 

Objective function: initially the weighted sum of squares of 
differences in frequencies; after a feasible design is 
found, the objective is changed to the weight. 

Design variables: flange thicknesses, lumped weights 

Boundary condition(s) at root: one analysis performed with pinned 

conditions, another analysis perfomed with fixed conditions | 

Frequency Constraints: windows on first, second and third 

collective flapping and also on first, second and third 

cyclic flapping I 

Autorotation constraint: yes | 

Stress constraint: no 

References: Third Semi-Annual Report pp.7-9, Thesis, pp. 50-54 

* 

— — — — — — —— — — — — — — — —— — — — — — — — j 

Case 13. Teetering rotor (similar to Cases 10-12, except that j 

collective and cyclic flapping and inplane and also '■ 

torsional modes considered) f 

j 

Rotating: yes 1 
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Objective function: initially the weighted sum of squares of 
differences in frequencies; after a feasible design is 
found, the objective is changed to the weight. 

Design variables: flange thicknesses, wall thicknesses of 

both sides of box cross-section, lumped weights, stiffness 
of torsional spring at root 

Boundary condition(s) at root: a) flapping — one analysis 

performed with pinned conditions, another analysis perfomed 
with fixed conditions; b) inplane — one analysis performed 
with pinned conditions, another analysis perfomed with fixed 
conditions; c) torsion — fixed conditions 

Frequency Constraints: windows on first, second and third 

collective and cyclic flapping; windows on first, second and 
third collective and cyclic inplane; and window on first 
torsional 

Autorotation constraint: yes 

Stress constraint: yes 

References: Third Semi-Annual Report pp. 12-13, Thesis, pp. 58-60 

Case 14. Teetering rotor (similar to Case 13, except that box- 
beam dimensions are fixed) 

Rotating: yes 

Objective function: weighted sum of squares of dirferences in 
frequencies 

Design variables: lumped weights, stiffness of torsional spring 
at root 

Boundary condition(s) at root: a) flapping — one analysis 

performed with pinned conditions, another anaxysis perfomed 
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with fixed conditions; b) inplane — one analysis performed 
with pinned conditions, another analysis perfomed with fixed 
conditions; c) torsion — fixed conditions 
Frequency Constraints; windows on first, second and third 

collective and cyclic flapping; windows on first, second and 
third collective and cyclic inplane; and window on first 
torsional 

Autorotation constraint: yes 
Stress constraint: yes 

References: Third Semi-Annual Report pp. 13-14, Thesis, pp. 60-61 

Case 15. Teetering rotor (similar to Case 14, except that 
stiffness of blade cross-section at root is a design 
variable) 

Rotating: yes 

Objective function: weighted sum of squares of differences in 
frequencies 

Design variables: lumped weights, stiffness of torsional spring 
at root, variable root-stiffness (but except at root, all 
other dimensions of the box cross-section are fixed) 

Boundary condition(s) at root: a) flapping — one analysis 

performed with pinned conditions, another analysis perfomed 
with fixed conditions; b) inplane — one analysis performed 
with pinned conditions, another analysis perfomed with fixed 
conditions; c) torsion — fixed conditions 
Frequency Constraints: windows on first, second and third 

collective and cyclic flapping; windows on first, second and 


third collective and cyclic inplane; and window on first 
torsional 

Autorotation constraint: yes 

Stress constraint: yes 

References: Third Semi-Annual Report pp.14. Thesis, pp. 61 

Case 16 . Teetering rotor (similar to Case 1% except that blade 
pretwist is included. 

Rotating: yes 

Pretwisted Blade: yes 

Objective function: initially the v/eighted sum of squares of 

differences in frequencies; after a feasible design is 
found, the objective is changed to the weight. 

Design variables: flange thicknesses, wall thicknesses of both 

sides of box cross-section, lumped weights, stiffness of 
torsional spring at root 

Boundary condition(s) at root: a) f lapping — one analysis 

performed with pinned conditions, another analysis performed 
with fixed conditions; b) inplane — one analysis performed 
with pinned conditions, another analysis performed with 
fixed conditions; c) torsion — fixed conditions 

Frequency Constraints: windows on first, second and third 

collective and cyclic flapping; windows on first, second and 
third collective and cyclic inplane; window on first 
torsionaT 

Autorotation constraint: yes 

Stress constraint: yes 

References: Thesis, p. 61, 64 
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Case 17. Teetering rotor (similar to Case 14, except that blade 
pretwist is included 
Rotating: yes 

Pretwisted Blade: yes 

Objective function: weighted sum of squares of differences in 

frequencies 

Design variables: lumped weights, stiffness of torsional spring 

at root 

Boundary condition(s) at root: a) flapping — one analysis 

performed with pinned conditions, another analysis performed 
with fixed conditions; b) inplane — one anlaysis performed 
with pinned conditions, another analysis performed with 
fixed conditions; c) torsion — fixed conditions 
Frequency constraints: windows on first, second and third 

collective and cyclic flapping; windows on first, second and 
third collective and cyclic inplane; window on first 
torsional 

Autorotation constraint: yes 

Stress constraint: yes 

References: Thesis, p. 61, 64 

Case 18. Articulated rotor 
Rotating: yes 
Pretwisted Blade: yes 

Objective function: initially the weighted sum of squares of 
differences in frequencies; after a feasible design is 
found, the objective is changed to the weight. 

Design variables: flange thicknesses, wall thicknesses of 

13 
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both sides of box cross-section, lumped weights, stiffness 
of torsional spring at root 

Boundary condition(s) at root: a) flapping — pinned; b) inplane 
— pinned at same radial location as in the case of 
flapping; c) torsion — fixed conditions 
Frequency Constraints: windows on first, second and third 

flapping; windows on first, second and third inplane; and 
window on first torsional 
Autorotation constraint: yes 
Stress constraint: yes 

References: Fifth Semi-annual Status Report pp. 8-9, 

Thesis, pp. 67-69 


Case 19. Articulated rotor (similar to Case 18, except that box- 
beam dimensions are fixed) 

Rotating: yes 
Pretwisted blade: yes 

Objective function: initially the weighted sum of squares of 
differences in frequencies; after a feasible design is 
found, the objective is changed to the weight. 

Design variables: lumped weights, stiffness of torsional spring 
at root 

Boundary condition(s) at root: a) flapping — pinned; b) inplane 
— pinned at same radial locacion as in the case of 
flapping; c) torsion — fixed conditions 
Frequency Constraints: windows on first, second and third 

flapping; windows on first, second and third inplane; and 
window on first torsional 
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Autorotation constraint: yes 
Stress constraint: yes 

References: Fifth Semi-Annual Status Report, pp. 9-10 
Thesis, pp. 67-70 

Case 20 . Articulated rotor (articulation at different stations for 
flapping and inplane motion) 

Rotating: yes 

Pretwisted Blade: yes 

Objective function: initially the weighted sum of squares of 

difference in frequencies; after a feasible design is found, the 
objective is changed to the weight. 

Design variables: lumped weights, stiffness of torsional spring at 

root. 

Boundary conditions (s) at root: 

a) Flapping — pinned; 

b) Inplane — pinned, but pin location is another few feet away 
from root 

Frequency constraints: windows on all first, second and third fre- 

quencies (flapping, inplane, torsion) 

Autorotation constraint: yes 

Stress constraint: yes 

References: Sixth Semi-Annual Status Report, pp. 3-4 



Case 21 . Articulated rotor (similar to case 20 except that box beam 
dimensions are also design variables) 

Rotating: yes 

Pretwisted Blade: yes 
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Objective function: initially the weighted sum of squares of differ- 

ences in frequencies; after a feasible design is found, the 
objective is changed to the weight. 

Design variables: flange thicknesses, wall thicknesses of both sides 

of box cross-section, lumped weights, stiffness of torsional 
spring at root. 

Boundary conditions (s) at root: 

a) Flapping — pinned; 

b) inplane — pinned, but pin location is another few feet awai 
from root 

I 

1 

Frequency constraints: windows on all first, second and third fre- 

quencies (flapping, inplane, torsion) 

Auto rotation constraint; yes 

Stress constraint: yes -* 

References: Sixth Semi-Annual Status Reports, p. 5 

3.2 FINDINGS RELATED TO OPTIMIZATION OF ROTOR BLADES 

' . t 

V/ 

The most important general finding of the project is that it • 

is possible to use an optimization routine such as CONMIN to 
tailor blade mass and stiffness distributions in an optimal 
manner. Furthermore, formulating the optimization problem in j 

i 

terms of frequency placement (that is, restricting the natural ! 

j 

frequencies of the blade to lie within narrow intervals located ] 

away from certain integer multiples of the rotor speed) has been \ 

j 

shown to be a useful approach for reducing vibrations. j 

, 5* < 

In addition to these general findings, the project estab- ^ 
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lished a number of specific results, knowledge of which would be 
useful to anyone intending to apply or extend the optimization 
approach developed during the project. A list of these results 
follows. 

1. In applying CONMIN to rotor-blade design, gradients of 
the objective and constraint functions should be calcu- 
lated by analytical formulas rather than by finite 
differences. However, finite differences serve as a 
useful check on the possibility of errors in the com- 
puter implementation of the analytical formulas. 
Reference: Second Semi-Annual Status Report, pp. 9-10; 
Thesis, pp. 23. 

2. Frequency constraints may be formulated directly in 
. terms of the frequency in Hz , rather than in terms of 

eigenvalues (i.e., the square of the circular frequen- 
cy). If eigenvalues are used, then scaling should be 
employed in the constraint equations to ensure well- 
behaved gradients for use in CONMIN, 

Reference: Second Semi-Annual Status Report, pp. 10-12; 
Thesis, pp. 25 

3. The following values of CONMIN parameters were ade- 
quate for most of the optimization studies: 

ITMAX » 40-80 
ITRM = 3 

DELFUN = 0.0001 (for cantilever beams) 

* 0.00001 (for rotor blades) 
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DABFUN = 0.0025 (for cantilever beams) 

= 0.00001 (for rotor blades) 

THETA = 1.0 
PHI = 5.0 

Reference: Second Semi-Annual Status Report, pp. 13-16; 
More efficient designs can be achieved if lumped 
weights are included as design variables (along with 
dimensions of the cross-section of the blade). 
Reference: Second Semi-Annual Status Report, pp. 18-21; 
Thesis, pp. 27. 

Because of the stiffening effect of the centrifugal 
forces in a rotating blade, frequency placement is much 
less dependent on stiffness and mass distributions than 
in a non-rotating blade. Thus, the rotational speed has 
a strong influence on what can be achieved in the 
optimization process. 

Reference: Second Semi-Annual Status Report, 
pp. 23-24; 

Use of ten finite-elements appears adequate to model a 
rotor-blade for optimization studies, although if many 
frequencies must be calculated, more elements must be 
used. Empirical rules which have been suggested are a) 
Use 4n degrees-of- freedom; and b) use n 4 degrees, where 
n is the number of frequencies to be found. 

Reference: Second Semi-Annual Status Report, pp. 25-29; 
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An accurate eigenvalue routine should be used in opti- 
mization studies, since errors in eigenvalue calcula- 
tions can appreciably affect the optimal design. 
Reference: Second Semi-Annual Status Report, pp. 29-31; 

The natural frequencies of a rotor blade are not great- 
ly affected by small changes in dimensions of the blade 
cross-section . 

Reference: Second Semi-Annual Status Report, pp. 31-32; 

For tight frequency-constraint windows, CONMIN is often 
unable to find a feasible design. In such cases, an 
objective function consisting of the weighted sum of 
squares of the differences in frequencies (actual fre- 
quency minus desired value) may be used initially. In 
the process of minimizing this objective function, 
CONMIN is often able to find a design which satisfies 
the frequency constraints. If this occurs, the objec- 
tive function may then be switched to the weight of the 
blade . 

Reference: Third Semi-Annual Status Report, pp. 5-9; 
Thesis, p. 44, 48, 50, 52, 58 

The natural frequencies of a blade which has already 
been built can be modified in a rational (rather than 
trial-and-error) manner by using CONMIN to specify 
where lumped mass should be added. It appears best to 
either raise all undesirable frequencies or lower all 
undesirable frequencies (a mixture of raised and 


* rW 
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lowered frequencies is much more difficult to attain). 
References: Third Semi-Annual Status Report, pp. 13-14; 
Thesis, pp. 60-61; Sixth Semi-Annual Status Report, 
p. 4. 

The forced response of a rotor blade can be adequately 
controlled through the approach of frequency placement. 
Reference: Fourth Semi-Annual Status Report, pp. 2-10; 
Thesis, pp. 71-95. 

Aerodynamic damping substantially reduces resonant 
peaks, but even in the presence of damping, frequency 
placement is a powerful driver of loads, and, as a 
result, frequency placement can be justifiably con- 
sidered an important part of blade optimization in the 
presence of damping. 

Reference: Fourth Semi-Annual Status Report, pp. 5-7; 
Thesis, pp. 75-82. 

Finite-element modelling errors caused by neglecting 
secondary structural items such as shear deformation, 
restraint of warping during twist, and filler stiffness 
are small. However, accurate filler properties, dimen- 
sions and locations are required in order to model the 
mass distribution properly. 

Reference: Fourth Semi-Annual Status Report, pp. 11-12; 

Since calculating eigenvalues is the major computa- 
tional burden in rotor-blade optimization, an efficient 
eigenvalue routine should be used. For example, determ- 
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inant search or subspace iteration can be used to 
calculate only the needed first few frequencies. 
Reference: Fifth Semi-Annual Status Report, pp. 2-5 
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The mass and stiffness distributions for helicopter 
rotor blades ar e to be tailored in such a way to give a 
predetermined placement cf blade natural frequencies. The 
optimlal design is pursued with respect of minimum weight, 
sufficient inertia, and resonable dynamic characteristics. 
Finite element technique will be used as a tool. Rotor 
types include hingeless, articulated, and teetering. 
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DESIGN OP HELICOPTER ROTOR BLADES FOR 


OPTIMUM DYNAMIC CHARACTERISTICS 


1. INTRODUCTION 

1.1 HELICOPTER DESIGN 

The design of helicopter rotor blades involves not only considera- 
tions of strength, survivability, fatigue, and cost, but also requires 
that blade natural frequencies be significantly separated from the 
fundamental aerodynamic forcing frequencies (e.g. Ref. 1). A proper 
placement of blade frequencies is a difficult task for several reasons. 
First, there are many forcing frequencies (at all integer-multiples of 
the rotor RPM) which occur at rather closely-spaced intervals. For 
example, J/rev and 6/rev are less than 20 % apart. Second, the rotor RPM 
may vary over a significant range through the flight envelope, thus 
reducing even further the area of acceptable natural frequencies. Third, 
the natural modes of the rotor blade are often coupled because of pit^h 
angle, blade twist, offset between the mass center and elastic axis, and 
large aerodynamic damping. These couplings complicate the calculation of 
natural frequencies. In fact, the dependence on pitch angle makes fre- 
quencies a function of loading condition, since loading affects collec- 
tive pitch. Fourth, the centrifugal stiffness often dominates the lower 
modes, making it difficult to alter frequencies by simple changes in 
stiffness or mass. 

In the early stages of the development of the helicopter, it was 
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believed that helicopter vibrations could be reduced (and even eliminau- 
ed) by the correct choice of structural coupling and mass stiffness 
distribution. However, it is easy to imagine how difficult it is to find 
just the proper paraaeters such that the desired natural frequencies can 
be obtained. The difficulties in placenent of natural frequencies have 
led. in many cases, to prellainary designs which ignore frequency place- 
ment. Then, after the structure is 'finalized' (either on paper or in a 
prototype blade), the frequencies are calculated ( or measured) and 
final adjustaents made. Reference [2] describes the development of the 
XH-17 helicopter in wfaloh a 300-lb weight was added to each blade in 
order to change the spanwlse and chordwlse mass distribution and therby 
move the first flapwise frequency away from 3 /rev. The authors were con- 
fident that similar adjustments to the mass distribution (and thus to 
the frequencies and modes of the blades) could greatly reduce rotor 
vibration on other rotors. An analytic study in Reference [3] predicts 
that chordwise mass distribution could also be used to lower overall 
helicopter vibrations. In particular, a forward shift of mass is shown 
to be useful because it places torsion in resonance with a particular 
harmonic. The torsion loads car. then be tuned to cancel undesirable 
blade loads. The study also shows, however, that such mass changes may 
have an adverse effect on stability; and thus stability and vibration 
must be studied together. Similar benefits of inertia pitch-flap coupl- 
ing are also observed in shaker tests In Reference [4]. 

These positive results, and other like them, were at least partial- 
ly responsible for the optimistic outlook so aptly presented in Ref [51. 
In that reference, six helicopter pioneers express their belief that 
helicopter vibrations can be reduced through proper blade and fuselage 
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design. This optimism of the 50 's was somewhat eroded In the 50 's and 
70’ s as the true complications of rotary-wing dynamics became bettor 
lenovn. Nevertheless, the belief is still held by most dynamlclsts that 
simple concepts (such as frequency placement) can so a long way toward 
improving rotor design. For example, in Reference [51 six helicopter 
pioneers (some of the authors of Ref. 4), reminisce on the early days of 
rotary wing and the recent advances in our understanding of helicopter, 
let, they still contend that much can be learned from simple principles. 

Presently, helicopter blades are not tailored to give a set of 
desired natural frequencies. Instead, blades are designed based on other 
considerations (including the desired aerodynamic characteristics and 
the cumulative experience of the designers) . Then, after the design is 
analyzed (either by computer program or by fabrication and testing), the 
designer checks for frequencies that are poorly placed. These are then 
adjusted by judicious application of lumped inertias at crucial spanwise 
locations. These after-the-fact alterations, however, can be detrimental 
to blade weight, blade cost, and the development time of the aircraft. 
Sometimes, the problems are unsolvable, and a helicopter is left with a 
noticeable resonance problem. 

The state-of-the-art in helicopter technology is now to the point, 
however, that it should be possible to correctly place rotor frequencies 
during preliminary design stages. There are several reasons for this. 
First, helicopter rotor blades for both main rotors and tail rotors are 
now being fabricated from composite, materials (Refs. 7 and 8). This 
implies that the designer can choose, with limited restrictions, the 
exact El distribution desired. Furthermore, the lightness of composite 
blades for the main rotor usually necessitates the addit.n of weight to 



give sufficient autcrotational Made inertia. Thus, there is a consider- 
able amount of flexibility as to how this weight may be distributed. 
Second, the methods of structural optimization and parameter identifica- 
tion are now refined to the point where they can be efficiently applied 
to the blade structure. Some elementary techniques have already been 
used for the design of rotor fuselages (Ref. 9) . It follows that the 
time is right for the use of structural optimization in helicopter blade 
design. Some work on this is already under development, and, although 
not published, some companies are already experiment with the optimum 
way to add 7 weight to an existing blade in order to improve vibrations. 
1.2 PREVIOUS WORK 

In this light, we would like to mention a few recent attempts at 
application of optimization techniques to rotor blade design. In Ref 
[10], an optimization procedure is applied in order to reduce blade 
loads consistent with aeroelastic stability. The procedure is not com- 
pletely automated, however, and the designer must make the design incre- 
ment at each iteration based on numerical sensitivity parameters. The 
biggest needs (as identified in this work) are the complete automation 
of the optimization and the formulation of realistic design constraints. 
In Reference [11], an optimization package is applied to an aeroelastic 
response program. The results mirror the earlier conclusions of Refer- 
ence 2-4. In particular, minimization of vibrations tends to drive some 
natural frequencies close to Integers in order to cancel loads. (The 
rotor becomes an Isolator.) Although this turns out to be a good mathe- 
matical solution, stability analyses in Reference [11], as in References 
[2-4] . show that the coalescence of frequencies to suppress vibration 
tends to introduce aeroolastic instabilities. Thus, flutter margins tend 
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to bacon the dominant constraints. Furthermore, minimization of loads 
at one flight condition may not at all minimize loads at others. 

Another investigation into vibration reduction by alteration of 
mass and stiffness distribution is given in Reference [12]. In that 
reference, a tip weight is used to change the mode shape. It is hypothe- 
sized that changing the mode shape such that it is orthogonal to the 
forcing function is a way to lower vibrations. However, ths conclusions 
are uncertain since the frequencies also are changed by this added 
weight (e.g. the second flap mode moves away from 5.06 to 5.19 per rev). 
One also notices that the loading distribution changes with flight 
condition so that modal shaping may help one condition but hurt others. 
Other related previous work is found in Reference [13]. That paper shows 
that design to minimum loads can reault in a disjoint solution. Fortu- 
nately, in helicopter problems we generally begin with an adequate (but 
not perfect) blade design. Thus, many questions such as this (i.e. 
disjoint in the design space) are automatically avoided. We already have 
a good first guess and merely wish to refine it. 

1.3 SC n PE 

In this paper we undertake a much less ambitious aim than the mini- 
mization of hub loads. Instead we look at the problem of using optimiza- 
tion techniques in order to place natural frequencies. Even within this 
reduced problem there are varying levels of complexity. For example, one 
could consider the retrofit problem: 

' Given a blade design find the amount and location of added 
masses required to move frequencies away from integer 
resonances . ' 

One could also consider th< basic design problem in which both stiffness 
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and mass distributions may be ohosen. In this paper, we treat this 
latter design problem in the lorm of three uncoupled problems: flap, in- 
plane, and torsion. 

The scope of this present work is not Just to find a mass and 
stiffness distribution to give desired frequencies. It is also to deter- 
mine meaningful constraints and objective functions that will result in 
realistic designs. In this area, several items are noteworthy. First, 
there is the airfoil envelop. Whatever the structural engineer designs 
must )ie within the airfcil cross-section. Second, there is mass balanc- 
ing. The center of mass of each section should be forward of the one- 
quarter chord. Third, there is the autorotational constraint. The blade 
must have sufficient mass moment-of-inertia to insure a safe autorota- 
tional capability. Fourth, there is strength. The blade must be strong 
enough to endure the centrifugal loads as well as the oscillatory bend- 
ing loads. This last criteria is the most elusive of the four. Designers 
know how to make a very soft section (hinge or flexure) which neverthe- 
less can withstand high centrifugal and bending loads. Such flexures 
generally do not fall within an airfoil envelope, however, and are 
placed near the root. Therefore, in the work to follow, we first obtain 
'optimum' designs and then check to see if the required El distribution 
has unrealistically soft spots. Similarly, we check the final designs 
for axial stresses. 

In summary, we work with simple (buc realistic) rotor-blade designs 
and simply experiment with constraints and objective functions in order 
to determine the feasibility of designing to a ierired set of frequen- 
cies. 
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1.4 OVERVIEW OP OPTIMAL STRUCTURAL OBSIGN 

Most approaches to optimal structural structural design may be 
classified into three categories. (For . eeent review articles see Refs. 
14 aud IS.) One such category is 'variational methods.' These generally 
rely on techniques from the mathematical theory of the calculus of 
variations, and, when appllc.'*>le, often provide useful physical insight 
into the nature of an optimal design. Unfortunately, only relatively 
simple problems can be solved by this approach, since the mathematics 
becomes Intractable when complex engineering structures are considered. 

A second category of structural optimization techniques consists of 
the application of mathematical programming methods together with the 
discretization of the structure by finite element techniques. This 
approach to optimization was founded in I960 (Ref. 16) with the hore 
that more complex structures could be analyzed than were possible when 
using the analytical tecnniques of the calculus of variations. Rowever, 
in the late 60 's it became apparent that mathematical programming mothod 
had limitations of their own, namely, unacceptably long computation 
times occurring when the number of design variables become large (over 
20-100, depending on the type of structure). Fortunately, several 
Improvements developed over the last few years appear to have 
significantly extended the capability of the mathematical programming 
approach, and, as a result, it is this approach we intend to draw upon 
for solution techniques in this research. 

A third eatagory . f structural optimization approaches is the 
’optimality criterion' approach in which an equation expressing some 
necessary condition of optimality is used as the basis for constructing 
an iterative (successive re-design) procedure. Originally developed 
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because of dissatisfaction with mathematical programming techniques, the 
optimality-criterion approach Initially relied on intuitive optimality 
criterion such as constant stress-ratio and uniform strain-energy 
density conditions. More recently, optimality criterion (and associated 
re-design equations) have been derived from the Kuhn-Tuelcer conditions 
(see, e.g. Ref. 17) for a constrained minimization problem. 

The optimality criterion approach seems- especially well-suited to 
problems wi*" s . a large number of design variables. Since our design 
problem will have a moderate number of variables and since deriving 
efficient re-design equations for our problem is not immediately 
straightforward, we initially prefer tSe mathematical programming ap- 
proach orer the optimal! tv-oriterun approach. 

A structural optimization computer program, called CONMIN (Ref. 
27), is available from NASA. It is this program that is used in our 
present work. CONMIN is based on the mat^ciiatical nonlinear programming 


method of feasible directions. 



2. 3ACE GROUND 


2.1 FORMULATION OF PROBLEM 

Because numeric all/— baaed optimization is best carried out with dis- 
crete variables » the finite element technique stands as the most logi- 
cal choices for the blade model. A recent research project (Ref. 18) has 
resulted in a finite-element computer program that, is iu-jaliy suited to 
the work here. The program allows for tapered, twisted finite elements 
in a rotating environment. The existing code can calculate natural fre- 
quencies, (with ard without aerodynamic terms) and force response. 

Another important aspect of the rc*or blade optimization problem is 
the selection of the optimality criteria and constraints to be impose/ . 
Cur design problem has certain features which are unusual compared to 
typical problems occurring in the structural optimization literature. 
There are basically three catagories of criteria. In the first class. one 
would minimize weight given constraints on the natural frequencies (i.e. 
frequency ’windows'). In this case, a constraint on rotary inertia is 
also implied since a rotor must have sufficient inertia to autorotate. 
The advantage of this approach is that it is directly related to the 
physical realities of design. The disadvantage, however, is that the 
first guess will probably not be feasible (that is will not have fre- 
quencies that fall in the 'windows'). This can be a stumbling block to 
convergence. A second type of criteria is one in which the objective ! 
to minimize the discrepancies between desired frequencies and actual 
frequencies. The constraint then becomes a window on autorotational 
inertia. Although this avoids unfeasible solutions, it does not direct- 
ly minimize weight (although weight is limited by the autorotational 
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constraint), An objective function oan be constructed that combined com- 
bined blade mass and frequency placement, but the relative weightings of 
the two components is not obvious. The third category of constraint is 
to minimize vibrations directly without regard to frequency placement. 
Although this appears on the surface to be the perfec* solution, there 
are problems. First, calculation of vibrations is an order-of-magnitude 
more difficult than the calculation of frequencies. Second, past efforts 
at this have resulted in strange designs, incompatible with standard 
helicopter practice. Third, there is still the problem of the weight- 
vibration trade-off. In this work, we intend to concentrate on the first 
two categories with some attention to the third. 

Another type of constraint involved in the problem is the limita- 
tion on structural properties. The blade planform. airfoil, and twist 
are chosen by the aerodynamicist on the basis of performance. The 
structural engineer must choose his design to “it in the aerodynamic 
envelope given. There are five structural parameters to be chosen: 1) 
flapping stiffness, 2)inplane stiffness, 3) torsional stiffness, 4) 
mass, and 5) torsional moment of inertia. In practice, these cannot be 
chosen completely independently. Figure 1 shows the envelope of a 
typical blade section. All stiffness is assumed to reside in a box-beam 
of dimension b x h with thicknesses t,dl,d2. This beam is placed as far 
forward as possible (to keep the elastic axis near the 1/4 chord). Mass 
properties are due to the box-beam, skin, honeycomb, and two lumped 
masses. The lumped mass in the tip is typical of rotor blades and is 
used to keep the mass center forward of the aerodynamic center. A second 
mass is included to allow Independent choice of mass and mass-moment. 


The constraints of this construction are clear and are listed on the 
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figure. 

In addition, there ’re minimum constraints on t,dl,d2 to hold 
centrifugal loads and to remain within manufacturable limits. For 
example a simple minimum constraint on area could come from the centri- 
fugal constraint (not considering bending stress). Thus, if ^ is the 
maximum stress and if f is a safety factor, then 
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Of course, when we enter the vibratory- response phase of the work, bend- 
ing stresses will be included. 

Our work will nevertheless Include flutter criteria in a simplified 
manner. First, we can choose frequency placement such that no coales- 
cence occurs between flap-lag, flap-torsion, or lag-torsion. Second, we 
can constrain the five parameters in Figure 1 such that the mass carter 
is always forward of the 1/ 4-chord, a common design practice to prevent 
torsion-flutter in rotor Hades. 

2.2 F ini t o-gl ement Mode 

Although tapered, twisted elements are within our capabilities, we 

introduce here a simple case which is also of value. The stiffnesses GJ, 

El , El are assumed to be constant along the length of the element, 
zz 77 
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The lumped mass weight is assumed to be evenly distributed on the two 
nodes. 

lat the deflection of an element in the y a. ' z directions at a 
distance x be denoted as w(x) and v(x), for which the displacement 
models are assumed to be polynomials of third degree. The expressions 
are given as 


w(x) .71 ( 2x3 - 3 ' x2 * ' 3 ) t 3 ' x2 - Zx 3 ) 

" ~2 ( “ 2/x 2 + / 2 x ) - -t ( x 3 - / x 2 ^ C2a) 

r 1 2 

v(x) = ( Zx 3 - 3/x 2 + f 3 } * -S- (3/x 2 - 2x 3 ) 

r 1 3 

♦ ^(x 3 -2/x 2 W 2 x) + j|(x 3 -,x 2 ) C2b) 


where v^, v 3> v^ and v g represent the bending degrees of freedom in the 
zx plane and u^, u ? . u^ and u ? represent the degrees of freedom in the 
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i) Tha strain energy due to bending deformation can be expressed as 
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11) The potential energy in tension from the centrifugal force field, 
which is equivalent to the negative of kinetic energy due to radial dis- 
placement, is given by 






(4) 


where T, tension force, is assumed to be constant along each ele- 
ment. 

ill) The kinetic energy due to inplane displacement is given by 


r • 1 ? t 

V ” 2 V ~ fr dx 


( 5 ) 


which is equivalent to u 


T 





- 15 - 


U 


1 


U- 

3 



u 


3 


U 


4 


U 


2 



Oegrses of Fraedom of an Element 


Meanwhile, the pretwiat angle <p (x) , and the torsional deformation 
Q (x) are assumed to be polynomials of first degree, and can be express- 
ed as 


♦W * +,(■ 1 - f > - ♦,<- 7 ! 
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where <p t , <p x represent the pretwist angle at node 1 and 2, and Uj,u 10 
represent the elastic torsional degree of '’reedom at each end. 
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iv) The torsional energy, due to elastic deformations and centrifugal 
terms, can be expressed as 
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v) The 'torsion-rotation' energy under the effect of rotation is 
given by 


(km 2 len| )(<*> + d } 2 dx ( 

2 2 

where taa^ , km^ are mass moment of Inertia which can be expressed as 

km 2 * //pz 2 dy dz - P I ( 

^2 * ff pyZ dy dz 3 p r :z ( 



Total displacement energy now can be used to form the stiffness matrix 
from 



(3) 


r 

a 



•sir 


where u is the rector of nodal displacements. In the order as u^, 
u t , u g , u,, u^, u^, u 9 ,Uj ,u 1Q ; [K] is the elemental stiffness matrix 
of order 10. 

ri) The mass matrix will be obtained by the kinetic energy of an 
element, which is given by 
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Written in matrix form, the kinetic enerwy can be expressed as 
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where [M] is the mass matrix 
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3. ILLUSTRATIVE EXAMPLES OP STRUCTURAL OPTIMIZATION 

3.1 CANTILEVER BEAM WITH GIVEN FREQUENCY 

Soma staple examples will be examined and discussed before the 

utilization of the program CONMIN. In each case the results will be com- 
pared with those obtained by previous researches, if it is available. 

The first limiting example is the problem of determining the opti- 
mal design of an elastic Cantilever beam, such that with a specific 
natural frequency, the weight of the structure attains the minimum value. 

We start with a uniform beam, modeled by ten elements, with a given 
length of 1C inches, B ■ 1.0 lb-in 2 , El - 10 lb-in 2 , density - 0.042 
lb/in 2 , and a specified first lowest natural frequency ■ 0.6489 rad/sec. 
We obtain the final stiffness profile shown in Figure 2. Figure 2 is the 
present result with ten elements. 

A related problem has also been treated by Olhoff [19]. He 
seeks the design of a cantilever beam that yields a maximum value of a 
particular higher natural frequency w R (i.e. of specified order, n) 
with the volume and length of the beam specified. His work is the dual 
problem of the example shown in Figure 2. Optimization with respect to 
the frequency under the constraint of volume is similar to the one of 
minimizing weight (or volume) under the constraint of specified natural 
frequency. Figure 3 gives the profile of the optimal cantilever for n * 
1 by Olhoff. One can see that the shapes in Figure [2] and [31 are very 
simlllar in that they give a nonlinear taper. 

3.2 CANTILEVER WITH TIP MASS 
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Another example problem is to minimize the weight of a cantilever 
carrying a mass at the tip, subject to the constraint that the fundamen- 
tal natural frequency must be greater than or equal to a specified 
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ngure 2. Area Moment of Inertia for Optimum Beam in 
Present work. 
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Figure 3. Area Moment of Inertia for Optimum Beam 
from Reference 11. 
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value. The problem was originally formulated by Turner [201, Kahn anQ 
Willmert [21] used an optimality criterion method to solve Turner's 
problem. In this example, four fir ice elements are used, with the areas 
of each as the design variables, as illustrated in Figure 4. The spec- 
ified natural frequency is 17.752 rad/sec. The other initial data are 


Modulus of elasticity 
Mass density 
Radius of gyration ) 
Radius of gyration (A^ ) 
Radius of gyration (A 3 ) 
Radius of gyration (A^ ) 
Concentrated mass 
Length of each element 


- 10.3 x 10 psi 

- 2.5 x 10 lb-s 2 /in 

- 2.0 in 
*1.5 in 
* 1 in 

- 0.5 in 

- 1 lb-s 2 /in 
■ 60 in 


where I * Area * (radius of gyr-uion)' 
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The results of the optimization are shown in Table 1. The feasible 
starting design is described by k^ «* 200, k 2 * 15C, Aj - 60 and A^ - 35. 


Table 1 

Comparison of Can til ever Beam wl oh Concentrated mass 




Ref. [12] 

Ref. [13] 

This Paper 

Iteration 

- 

23 

10 

A 1 

(in 2 > 

136.81 

136.63 

134.60 

A 2 

(in 2 ) 

113.73 

118.70 

116. 

A 3 

(in„) 

L 

83.59 

83.58 

82. 

A 4 

(in.,) 

34.43 

34.61 

34.89 

Weight (lb) 

2243.0 

2242.9 

2214.4 

can be 

seen that 

excellent results 

* 

have been 

obtained using 


present COKMIN optimization program. 
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4. NUPlErtlCAL EXPERIMENTS 

4.1 DESIGN VARIABLE 

Despite the strong documentation and intensive development that has 
gone into optimization programs, it is always advisable to do some ex- 
perimentation with these programs for the particuiai class of problems 
to which they are to be applied. This has been done in detail for the 
representative box beam shown in Figure 5, The parameters ('see Fig 1) 
for this car. are: h » 2.5 in, s *0.1 in, s = 0.1 in. b = 4 in, t 
variable. This bean has been analyzed for various values of the CONMIN 
parameters and for various combinations of constraints. The first 
studies are performed for vertical vibrations, a nonrotating beam, no 
lumped mass, and with two frequency constraints. Each of the CONMIN 
options is then exercised, and several conclusions drawn. 

FiiSt, we find that the U3e of analytic gradients ithe derivative 
of objective function and constraints in closed form) is greatly to be 
desired. For the particular case in Figure 5, area, weight, and moment 
of inertia can be expressed in terms of the single variables, t 

A i - 0.5C + .60 t i (12a) 

V i - 1.20 ♦ 18.24 t L , (12b) 

- 25/90 + t i *(285 - 228 ^ + 60.8 t i 2 )'*124 (12c) 

Therefore, analytic derivatives are straightforward . Where analytic 
gradients are not available, however, we find that finite difference 
gradients still work albeit at a higher computationax cost. Second, wa 









find that the optimization is best behaved when frequency constraints 
are provided in Hz. A constraint (if not scaled) on eigenvalues (w ) is 
more difficult for the program when default values are used. Third, we 
find that initial designs outside of the desired constraints 
(infeasible) sometimes can lead to convergence. Since this is not always 
the case, however, alternative strategies are necessary. Fourth, we 
find that the default values for the CONMIN program worked reasonably 
well (although they are not always the most efficient values) . An 
example is the number of iterations. Sometimes 40 iterations were 
required for convergence, although the default value is 10. 

In terms of various modes of application, we also have cone to 
several conclusions. First, we could find optimun designs no matter how 
tightly we closed the windows on frequency (i.e. the frequency const- 
raints). Thus, we are able to essentially 'zero' an objective function 
based on frequencies (for vertical vibrations alone). Second, we can 
handle a large number of simultaneous frequency constraints. (We have 
successfully gone from 2 up to 5 constraints . ) The optimization also re- 
mains well-behaved when we add blade rotation, lumped mass, and the auto- 
rotational constraint. 

4.2 CONVERGENCE 

We have studied the convergence of the ^inal design as a function 
of the number of elements used in the finite-element frequency calcula- 
tion. Tc study how ths optimal design changes as the number of elements 
increases, a cantilever beam with 'n' elements and with lumped weights 
added at the nodes but otherwise similar to ths beam in Figure 1 is 
considered. The density, and the constraints on the natural fi equency, 
lurped weights, and moments of inertia are 



I 
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P - 0.05 lbf/in 3 

1.0 < f x < 1.3 (Hz) 

10.0 <. n - 11,2 (Hz) 
o.o < w 1 < ioo.o (ibfj 

3.83073 < < 5.2083 (in 4 ) 


(13a) 

(13b) 

(13c) 

(13d) 


and the initial value of OBJ (the total weight) is 30.2249 Ibf. 

Results of the study are shown in figures [6-llj. In all cases, 
the active frequency constraints were found to be 

f t - 1.3 (Hz) 
f 2 - 11.7 (Hz) 





Figure [81 demonstrates, as one would expect, that the optimum 

weight does in fact decrease as more elements are added to the mesh. The 
change in optima weight is quite snail (note that the scale of the 

vertical axis begins at 20.0). 

Figure [71 and [8] show the variation of the lumped weight and the 
moment of inertia (of the cross-sectional area) at the free end versus 
the total number of elements in the mesh. It appears that these quanti- 
ties do not converge. The result cam be explained, however, by referring 
to Figure [91, in which the upper curve represents the total weight at 
the free end. ( m(n) is the non-structural, or, lumped weight i(n) is 
the structural weight associated with the mass distributee* through' it 


element 'n'). it can be seen from the figure that the total weight 
appears to converge smoothly as the mesh is refined. The explanation for 
the apparent non-convergence shown in Figures [7] and [81 and for the 
convergence shown in the top curve of Figure [9] is that the 'structural 
weight' at the free end of the cantilever is not really structural, 
since there is no portion of the beam beyond the free end which needs to 
be supported. Thus, the optimization routine is indifferent to whether 
structural or non-structural weight is present at the free end — the 
only thing that counts is the total weight at that end. 

Figure [91 also shows the var'^tion of the lumped weight slightly 
beyond the middle of the beam. (All optimal designs have non-zero lumped 
weights there and at the free end of the beam.) The weight can be 
decrease smoothly as the mesh is refined, although no asymptote apoears 
to be present, Ths explanation for this behavior is that, as the mesh is 
refined, the weight in the middle Is being place more efficiently - and 
thus less is needed. 

The various sketches in Figure 10a show the distribution of mass and 

stiffness along the beam for increasing numbers of elements. It is in- 
teresting to observe that the optimization routine finds It most effici- 
ent to meet the constraints on frequency by varying the lumped weight 
rather than by varying the stiffness (moment of inertia) , since this 
latter quantity is at its lower bound every re except near the end of 

the beam. Another interesting aspect of Figure 10 is the manner in which 
the lumped mass at the center alternates between: 1) being all on one 

element, and 2) being split between two elements. This phenomenon is a 
result of the fact that the minimum weight structure would have all the 
mass at a single point (node or antinode). When this single point lies 


OBJ (LBF) 
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TOTAL NUMBER OF ELEMENTS IN MESH 

Figure S Convergence of Weight of Objective Function 
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near a structural node the mass is placed there. On the other hand, when 
the mesh causes the point to be between nodes, the mass is accordingly 
divided between the two closest nodes. 

As was pointed out previously in reference to Figures [73 and [8], 
the optimization algorithm appears to treat the structural and non- 
structural mass at the end of the beam as interchangable. To test this 
hypothesis further, the optimal design problem statement was altered 
slightly by decreasing the upper bound constraint on the moment of iner- 
tia from 5.2083 to 2.0. The resulting optimum design is shown in Figure 
[10b], and should be compared with the design (for n = 10) 3hown in 
Figure [10aJ,, Mote that the constraint on the moment of inertia for 
element 10 is not active in the optimal design of Figure [10b] (the 
constraint was active during the CONMIN iterations leading to this 
optimal design). Thus, the effect of the constraint is to lead the 
optimization algorithm along a different path than that followed when 
the constraint value was 5.2083. The design found, however, has about 
the same total weight at the free end (* 9.9705 lbf) as the previous 
ten-element optimum (* 9.9222 lbf). Thi3 result confirms the hypothesis 
that CONMIN increases the moment of inertia at the free end only as a 
means of increasing the mass there. Once that option is closed (that is, 
the upper bound constraint is reduced to a value of 2.0), CONMIN simply 
increases the l.cped weight at the beam tip. This finding suggests 
that, in future optimization studies, a tight constraint be imposed on 
the moment of inertia at the free end, since little structural capa- 
bility is needed there, and necessary end mass can be adequately repre- 


sented by the lumped weight design variables. 
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5. PRELIMINARY CALCULATION FOR ROTORS 
5.1 WIND TURBINE BLADE 

The first example In this section Is the optimization of a wind 
turbine rotor blade at 30 rpm. Initial data is taken from Ref. [22], A 
ten-element model is used. Only the flapping is considered. The area 
moment of inertia. I, and the lumped weight of each element are taken as 

the design variables (see Fig 1 for blade area cross-section). Young's 

6 2 3 

Modulus, E - 0.2 x 10 lb-in , and density * 0.0334 lb/in are assumed 

to be constant. Blade radius, R * 750 inches. Table 2 shows the profile 

of moment of inertia and the distribution of added weight for the 

initial and final configurations. The final profile of the area moment 

of inertia along the blade is similar as the one in the previous 

example. The optimization procedure has removed material from the 

inboard sections and placed it more outboard. The lumped mass is 

concentrated at the tip of the blade as might be expected. We also note 

that most of the originally-postulated lumped mass is removed so that 

only the mass inherent in the stiffness elements or necessary for the 

autorotational constraint is maintained. (Although wind turbine have no 

autorotational constraint, a certain moment of inertia is still useful 

to smooth out wind vibration. 

An important aspect of the optimization problem is the existence 
(or lack of it) of a feasible solution. A 'feasible solution' is defined 
as any set of design variables that satisfy the constraints (whether or 
not that particular solution is an optimum). Jt is possible that, if the 
problem is poorly formulated, no feasible solution exists. What is 


more often the case, however, is that there are feasible solutions but 
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that the optimization scheme may not be able to find them. Thus, it is 
advantageous to have a feasible Initial guess so that one is assured 
that at least a local optimum is possible. 

For example. Table 2 illustrates that the first guess is feasible 
(Wj> 2.82 no. /rev). Here we found that CONMIN was able to move from this 
first guess through the space of feasible solutions. In other cases, 
however, when the first guess is not feasible we have found that CONMIN 
is not able to reach a solution. In such cases, one must add or remove 
some weight (or add or remove El) from the first guess to move from ini o 
the feasible space; or, alternatively, we must begin with frequency- 
placement as the objective function and then switch to weight when the 
frequencies are within tolerance. 

For example. Table 3 represents data for the s— e wind turbine as 
in Table 2, but the constraint on the first natural frequency has been 
lowered to remove it from the dangerous 3/rev range. This implies that 
the first guess in Table 2 is no longer feasible. In order to overcome 
this, a lumped mass is added to station 9 (225.4 vs 49.50). This lowers 
Wj^ below 2.621 rev but also lowers w^ to 8,25/rev. This could be alle- 
viated in one of two ways: 1) move the mass to the node of the second 
mode, or 2) simply widen the w^ window. We have done the latter. It is 
interesting that the added weight is ultimately rearranged to other 
places and other weight removed such that the new design Is no heavier 
than the optimum in Table 2. Furthermore, is raised to 8.57 so that 
the 'widened window' had nc effect on tne solution. 





9 




S' . - ?<* 



Table 2. Optimization of Wind Turbine Blade 
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3.2 HELICOPTER BLADE 


The design and analysis of a representative helicopter olade is 
discussed. Similar to Section 3, only flapping is considered and a te. , 1 - 
element model is used. The initial configuration is modeled after the 

rotor in Reference 23. Density is constant along the blade and equal to 

-3 3 8 2 

0.17 x 10 slugs/ in . Young's Modulus is equal to 0.49 x 10 lb-in at 

g 

the root and is equal to 0.585 x 10 lb-in elsewhere. Blade radius is 
equal to 193 inches. Results are given in Table 4 and 5. In Table 4, w^ 

is in the desired range but w 2 is too small. Furthermore, the autorota- 
tional inertia is larger than necessary. In this case, the C0NMIN pro- 
gram is able to remove mass and stiffness in such a way to raise w 2 and 
lower w^. The minimum bending inertia set as a constraint (0.4) is 
reached at every point except the root. The root remained high to keep 
w^ > 1 05, The new blade is one-third the original mass. In Taole 5, a 
stiffer initial design is used and the frequency is forced to be 
very high. In this case, the program C0NMIN would 'like' tr decrease El 
and m, but any removal of material could lower w^ beyond its lower 
bound of 1.24/rev. To counter this, the optimization oheme adds El near 
the root (to maintain w^ > 1. 24/rev). Furthermore, the lumped mass 
necessary to maintain autorctational constraint is moved sUght? y ■'n- 
board to have less effect on (keep it high) but more effect on w^ 

(keep it low) . This example illustrates the physical soundness cf this 
optimization scheme. It d>'es the same things that a designer would do 
(given appropriate constraints) but in a more systematic manner. Thus, 
with proper constraints, optimization can prove a valuable tool for 


frequency placement 



®p*nr 





INITIAL 141.76 8.883 1.25 10.55 2.76 8.88 
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We have also examined the designs in Tables 3-5 with respect to 
axial stress due to centrifugal loads. In each case, the maximum stress- 
es after optimization are equal to or only slightly higher than the 
original stresses. 
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6. TEETERING ROTORS 


6.1 Definition 

In this section we attempt our firsc optimization of a realistic 
cross-section (see Figure 11), one which is little different from our 
first design pattern, (Figure 1). Therefore, it is sufficiently general 
that both the bending and torsional stiffnesses of some currently exist- 
ing blades can be matched. Using this generic cross-section and starting 
from an actural rotor blade design. (Ref. [24]), we have studied the 
possibility of moving natural frequencies away from resonances while 
simultaneously satisfying constraints on the following: stress, the size 
of lumped weights to be added, capability for autorotation, and thick- 
ness of the main structural member (the box beam). Because a teetering 
blade was considered, cyclic and collective modes of vibration were 
calculated independently by a change in the boundary condition at the 
blade root. In the initial phase of the 3tudy, we considered collective 
flapping modes first, then cyclic flapping, and finally combined 
collective and cyclic flapping. The results of these studies were 
favorable (i.e. we were able to change the frequencies in the desired 
manner and still satisfy the constraints). Building on these results, in 
the second phase of this section, we consider a more challenging problem 
which involves combined modes of collective flapping, cyclic flapping, 
collective inplane, cyclic inplane and torsional vibrations. 

In this section, the primary design variables are: 1) the wall 
thickness of the box beam and 2) lumped weights, that can be added at 
specified stations along the beam. In the final problem studied, the 
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wall thicknesses are taken as fixed, and only the lumped weights are 
allowed to vary. This situation corresponds to that encountered in 
practice when a blade has actually been designed and manufactured, but 
then found to have poorly placed frequencies — thus lumped weights are 
added at various positions along the beam to change the frequencies. We 
found that our optimization routine was able to handle this problem 
adequately, although the total weight of the beam could not be used as 
the objective function, as had been done previously. Instead, a 
'frequency placement' objective function was used. 

6.2 Flapping Frequencies 

The first set of optimization problems in the section is concerned 
with flapping response only. The starting design for the optimization 
procedure in each of the three cases studied is a typical metal-bladed 
teetering rotor with a diameter of approximately 24 feet. Ten finite 
elements are used to model the rotor; their lengths are given in Table 
6. The objective of the optimization is to minimize the total weight of 
the blade. The design variables are the wall thickness, t^, of the 
finite element representation of the box beam (the structural member in 
the rotor - see Fig. 1 ) and the lumped weight w^ associated with each 
finite element. The lumped weight is the sum of two components, a fixed 
component (representing the weight of the leading and trailing edge 
strips, honeycomb, skin and nose weight,) and a variable component 
(representing additional ncm-structural ma.ss which may be added at 
various positions along the length cf the rotor) to modify the dynamic 
behavior in a desired manner. The side conditions on the element thick- 


nesses are, in units of inches. 


0.00044 < t jL < 0.730 

The side conditions on the lumped weights consist of a lower bound only, 
which represents the fixed component of weight for each element and is 
given in Table 6 under the heading *w j n '. Mote that the element thick- 
nesses t^ are not given in Table 6; instead, the area moment of inertia, 
X, is represented. Using the dimensions given in Fig. 12, we can show 
that I related to the thickness by the equation 

I - 2.593t 3 - 7.780t 2 + 7.780t + 0.506" 

The moment of Inertia is given, rather than the thickness, to facilitate 
comparison with I the portion of the moment of inertia which is con- 
tributed by those parts of the cross-section other than the box beam. 

(Thus remains fixed as t. is varied.) Table 6 also contains the 
o l 

values of the box weight, which are calculated by multiplying the weight 
density of the box beam material by the cross-sectional area of the box. 
Thus, the box weight is not an independent design variable, but depends 
on the thickness t^. The box weight is included in the table to facili- 
tate comparison with the distribution of lumped weight. The constraints 
for the optimization are both the autorotation constraint, 

sum <w i ')r jL 2 -> 0.5567xl0 7 lb-in 2 

(where w^ is the total weight of element i, and r^ is the distance from 
the root to the center of the i-th finite element). The frequency con- 
straints will be described in subsequent sections of this thesis.Some 
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Figi.'-e 12 Dimensions snd Definition of Design Variable for 
Box Beaa Cross-Section 
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additional data which complete the problem description are the values of 
the elastic modulus, 0.105x10 lb/in 4 , the radius, 288.8 in, the 
rotational speed, 324 rpm, and the mass density of the box beam 
material, 0.000262 uiugs/in^. 

6.2.1 Collective Modes 

The initial problem to be considered is the optimization of the 
blade with respect to collective flapping modes only. Because we are 
studying a teetering rotor, the collective mode of flapping may be 
analyzed by imposing a fixed boundary condition at the root cf the 
rotor. The imposed frequency constraints are 


1.15 

< 

Pi 

< 

1.50 

( 16 a) 

3.40 


P2 

< 

3.60 

( 16 b ) 

6.40 

< 

p3 

< 

6.60, 

( 16 c) 


in which pi, p2 and p3 are the first three collective flapping mode 
frequencies non— dimens ionali zed by dividing by the rotor speed. 

The starting design for the optimization algorithm is given in 
Table 6 under the beading 'initial'. This initial design was chosen to 
correspond closely with an actual rotor blade; thus.it is not surprising 
to find that the design is infeasible with respect to the frequency 
constraints we have imposed. The optimization algorithm used in this 
study, CONMIN, supposedly permits an infeasible starting point and at- 
tempts to proceed from this starting point to a feasible point. However, 
for our design problems, this feature of C0NMIN failed to produce a 
feasible design after many iterations. As an alternative approach, we 


•» 


s 





formulated a preliminary optimization problem in which the previous 
objective function (weight) was replaced by a 'frequency-placement' 
objective: 

obj fp - (3.50 - p2) 2 + (6. 50 - p3) 2 ( 17 ) 

The numbers 3.50 and 6.50 are the average of the bounds of the frequency 
constraint inequalities which are violated by th~ initial design. The 
remainder of the optimization problem is the same as the original pro- 
blem, except that the constraints on those frequencies which appear in 
the frequency placement objective are omitted. CONMIN was applied to 
this preliminary problem. In the process of minimizing the preliminary 
objective, CONMIN was able to drive the frequencies sufficiently close 
to their bounds that a feasible design (with respect to the original 
problem) was obtained. At this point, the original objective function 
was reinstated and CONMIN applied once again. 

Table 6 gives the optimized design obtained by this two-stage 
optimization procedure, with the corresponding frequencies and the total 
weight. From the point of view of helicopter vibrations, the initial 
design of this blade is acceptable, since (except for the third mode) the 
number/rev is far away fx m even integer values. The third mode is, 
however, near 6.0/ rev. The frequency of the second mode does not satisfy 
the inequality constraints, but is not near an even integer multiple. 
Note that the final design moves the third frequency to 6.53, wftile 
keeping the other frequencies within the constraints. At the same time, 
the weight of the blade drops from 344.5 lb to 265.6 lb. 


Table 6 Initial and Final dealgn for collective flapping *K>de* (tecterliiR rotor) 
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6.2.2 Cyclic Modes 


The next problem to be studied is the optimization of the blade 
with respect to cyclic modes of flapping. The non-dimensim«lized fre- 
quency constraints are now 


0.90 

< 

pl 

< 

1.05 

2.40 

< 

P2 

< 

2.60 

4.40 

< 

p3 

< 

4.60, 


in which pi, p2 .rad p3 are the first three cyclic non-dimens tonali zed 
flapping-mode frequencies. For cyclic flapping modes, the boundary con- 
dition at the root corresponds to a pinned support. The frequency- 
placement objective waa agein chosen by noting which frequency con- 
straints were violated by the initial design. Noting the initial fre- 
quency values given in Table 7, we define 


obj fp - (2.50 - p2)‘ ♦ (4.50 - P 3) 2 

Table 7 gives the optimal design found by the two-stage optimization 
procedure with the corresponding frequencies. The final weight of the 
blade is shown to drop from 344.5 lb to 295.2 lb. 

6.2.3 Combined Collective and Cyclic Modes 

Next we consider the optimization of the beam wit’ .‘•spect to 
combined collective and cyclic modes. Thus, in each itero.'io; an analy- 
sis must be performed to find the frecuencies corresponding ‘.o a fixed 
boundary condition; and then another analysis must be pe formed t^ c ir.d 



Table > Initial and Final design for cyclic falpplng nodes (teetering rotor) 
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the frequencies for a pinned boundary condition. All other aspects of 
the design problem remain the same as before. The constraints on the 
collective flapping modes are 


1.15 

< 

Pi 

< 

1.60 

(17a) 

3.40 

< 

p2 

< 

3.60 

(l?b) 

6.30 

< 

P3 

< 

6.60 

(17c) 

constraints 

on 

the 

cyclic nodes are 



0.90 

< 

Pi 

< 

1.10 

(18a) 

2.40 

< 

p2 

< 

2,60 

(13b) 

4.40 

< 

p3 

< 

4.60 

< 13c) 


For this problem, weighting factors are introduced into the frequency- 
placement objective. 


obj fp - [2(3.50 - p2) 2 ♦ (6.45 - p3) 2 ] . + 

J r K collective 

[2(2.50 - p2) 2 ♦ (4.50 - p2 > 2 i 

cyclic 


(13) 


The results of the optimization are given in Tables 8. The cyclic 
inodes of the initial design are well-placed in the sense that they are 
not near odd integers/rev, but the third collective mode is near 6.0/ rev 
(the same as in the first example). Note that the final design moves the 


ToHe 8 Initial and final design for collective and cyclic # 1 np|>* **rt ■mice 
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third collective frequency to 6.33 while keeping the other frequencies 
in the 'safe' range. However, the weight of the blade only decreased 
from 344.5 to 33S.5 lb, in contrast to the previous example in which the 
weight decreased to 265.6. The difference is caused by the larger number 
cf frequency constraints in this example compared to the previrus 
example. 

6.3 SIMULTANEOUS FLAPPING INPLANE AND TORSION 

The next problem to be studied is that of optimizing a teetering 
rotor blade subject to the following simultaneous constraints on fre- 
quencies (non-dimensionalized by dividing by the rotor rotational 
speed) : 

1 collective flapping modes. 


0.5 

< 

Pi 

< 

1.5 

2.3 

< 

P2 

< 

3.7 

4.3 

< 

P3 

< 

5.7 

cyclic flapping modes 

i. 




0.5 

< 

Pi 

< 

1.5 

1.3 

< 

P2 

< 

2.7 

3,3 

< 

P3 

< 

4.7 

collective inplane modes. 



0.0 

< 

Pi 

< 

1.0 

4.3 

< 

p2 

* 

\ 

5.7 

14.3 

< 

p3 

< 

15.7 


•I 


m*r 



4. cyclic inp lane nodes end 


0.1 

< 

pl 

< 

1.5 

(23a) 

5.3 

< 

P2 

< 

6.7 

(23b) 

17.3 

< 

P3 

< 

18.7 

(23c.) 


5. torsional node 

3.3 < pi < 3.7 (24) 

(The first loner bound for cyclic inplane nodes, 0.10, was later repalc- 
ed by 1.0 in the problea formulations of the following sections). Be- 
cause we are considering a teetering blade, the collective-flapping and 
cyclic— inplane nodes can be node led by clanped boundary conditions at 
the root, while the cyclic-flapping and collective- inp lane nodes can be 
modeled by pinned boundary conditions. 

The elastic modulus, blade length, speed of rotation, and density 
of the box-bean material are unchanged from the values used before. In 
addition to frequency and autorotation constraints, the axial stresa is 
constrained to be less than 20,000 psi. The value of the bound in the 
autorotational constraint has been changed slightly to 0.5429xl0 2 ib- 
in 2 . 

In the problem described in section 3., the thicknesses, t., d. ^ 

and both the vertical and horizontal walls of the box beam (see 
Pig. 2) are allowed to vary - that is, are also design variables, with 
the following side constraints (in units of inches); 
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0.01 < t. < 1.0 


0.01 < d u < 2.32 


(25a, b,c) 


0.01 < d 2i < 2.32, i - 1,2 10. 


The initial value* of these variables are given in Table 9. In the 
problem* discussed in section 3.2 and 3.3, the box bean dimension* t^, 
dj^ and d 2 ^ are fixed at these initial values. 

Tables 9 also gives data defining both fixed and initial stiffness 

and inertia values of the blade. In the tables, 'I ' and 'I ' repre- 

ox oy 

sent the portions of the flapping and inplane area moments of inertia of 
the blade section which are independent of the design variables. '1^' 
and '1^ ' are the area moments of inertia of the box beam - thus func- 
tions of the (initial) values of the design variables, w^, d^, d, , and 
t.. and -re ths rutary inertias of the box beam with respect to 

flapping and inplane and are calculated simply by multiplying the mass 
density of the box-beam material by the area moment of inertias. M ^ and 
M 0 y are the contributions to the rotary inertia of the section which are 
independent of the design variables. Rote that since these contributions 
come from items with different densities, a single uniform value of 

density cannot be defined for the M terms. 

o 

Other initial inertia and stiffness properties are also defined in 
Table 9. As before, the lumped weights associated with the finite ele- 
ments are taken as design variables; but, to permit greater latitude in 
placing torcional frequencies and to match more closely the behavior of 
a true helicopter blade, a torsional spring is introduced at the blade 
root; and it stiffness is taken as a design variable. The side 
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constraints on the lumped weights are given in Table 9 under the heading 
'w . the side constraint on the torsional stiffness consists of the 
requirement that the stiffness be non-negative. 

Since a torsional mode is involved, special treatment is given to 
GJ, the torsional rigidity, which is a function of all variables includ- 
ing t, dj, d 2 , the area of the trailing edge, and the lumped mass. The 
procedure for calculating GJ is described in Appendix A. 

Table 9 also contains the contribution of the lumped weight to the 

2 

rotary inertia, which equals the lumped mass (w^/g) times (b/2) . This 
expression has been chosen to match the behavior of the true blade. It 
is assumed that the lumped weights cf the first two elements contribute 
nothing to the rotary inertia. 

6.3.1 Variable Box Dimension 

As was done with the optimization involving flapping only, a two- 
step optimization procr dure is used, which involves a frequency- 
placement objective followed by a weight objective. The frequency- 
placement objective has the general form 

Obi'* - .u. [*f. <p„. - p.) 2 J (25) 

in which th um is taken over those f-equeneies which are to be changed 
from their ...itial values to the desired values p*.. Values of p*. and 
the weighting factors w^ are given in Table 10 with the results of the 
optimization. Examination of the values of lumped weight given in Table 
10 slows that weight is concentrated at the tip of the blade because of 
the autorotational constraint. - ae stress at the first element is close 

to the stress constraint value of 20.000 psi. Thi total weight changes 



Table 10 Initial and final daaign for flapping, Inplane, toralonal node a 
of Teetering rotor blade with variable boa be*'a dlaenaion 
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insignificantly, from 345 to 339 lb. However, the root spring shows a 
significant effect on the torsional mode. The spring constant changes 
from 6.5xl0 5 to 4.10x10^ in-lb/radian. 

The placement of the frequencies is shown in Table 9, also. Note 
that the third collective mode of flapping, which was near 6.0/rev, 
moves to 5.67/rev; and the second cyclic mode of inplane, which was near 
7/rev, moves to 6.7/rev. Finally, the torsion mode, which wa3 3.87/rev, 
moves to 3. 4/ rev. 

6.3.2 Fixed Box Dimensions 

We next consider the oroblem of modifying a blade which has already 
been constructed, but which has been subsequently found to have inappro- 
priate natural frequencies. Since the blade is already built, the only 
way its dynamic behavior can be modified is through the addition of 
lumped mass and also through changing the root spring. Thus, in contrast 
to the problem of section 3.1, here the box-beam weight, the flapping 
and inplane-bending inertia, and the torsional rigidity are constant. 
All starting values and fixed parameters are the same for this problem 
as in the previous section. Finding a starting design which satisfies 
all of the frequency constraints is a difficult task for this problem, 
and thus the frequency-placement objective is the only objective 
function used; the second phauie (with weight as the objective function) 
is never reached. Values of the weighting factors and frequency bounds 
which appear in the objective function are given in Table 11. 

The final values of the lumped weights and axial stresses are shown 
in same Table 11, with frequencies, root-spring stiffness, and total 
weight. The most significant frequencies are the third collective flap- 
ping mode and the second cyclic inplane mode, which are seen to move far 



away from the undesirable integer/rev values. However, the first cyclic 
inplane frequency is near 1/rev. 

6.3.3 Variable Root Bending Stiffness 

The problem formulated in the previous section presented computa- 
tional difficulties in that it was found difficult to find a design 
which satisfied all the frequency constraints simultaneously when only 
lumped mass and the root spring were used as the design variables. This 
difficulty may be caused by the dominating influence of the root bending 

inertias. Thus, it seems reasonable to include the values of 'I ' at 

oy 

the root as one of the design variables. The frequency-placement 
objective is used throughout the optimization (the weight is not used as 
the objective), and all starting data and fixed parameters are given the 
same values as in Section 3.2. Values of the weighting factors and 
frequency bounds which appear in the objective function are given in 
Table 12 with the results of the optimization. 

The results of Tables 12 differ from those of Tables 10 and 11. The 
lumped weight changes at the first, fourth, and fifth elements. The 

stiffness of the root spring moves from 6.5x10* to 4.32x10* in- 

lb/radian. However, the most significant effect is the change of I 

(the bending moment of inertia at the root from 48.11 to 290.29 in 4 . 

Table 12 shows that all frequencies are placed in tne safe range. 

6.4 Effect of Pretwist 

In this section, we will do the optimization of a beam which is 
pretwisted. (Pretwisted blade Implies that the mot' .*s of flapping and 
inplane are ccupled). The procedure of analysis will be the same as in 
section 6.3. Data are identical to those of section 6.3 save that the 


Table 11 lnltinl ami (Inal dealgn (nr flapping, Inpla ie and toralon.il aodea 
o i Teetering rotor blade wltli boa bene ulr-nelon fixed 
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pretwisted angle Is new included. 

6 . 4.1 Variable Box Dimension 

A two-step optimization procedure is still used. The box beam 
dimensions and lumped weights as (well as the root spring) are consider- 
ed as design variables. The results are shown in Table 13. In a similar 
fashion as results in section 6.3. the weight is concentrated at the tip 
of the blade because of the autorotatlonal constraint. Root stress Is 
close to the stress constraint which is 20.000 psi. The total weight 

does not decrease. Instead, it increases from '.b to 3S4 lb. The 

5 * 

spring constant changes from 6.5x10 to 4.2x10 ir-lb/rad. The second 
collective mode of flapping moves from 5.89 to 5.67 no/ rev while second 
cyclic mode of inplane moves from 7.08 to 6.54 no/rev. 

6.4.2. Fixed Box Dimensions 

For a blade of existing constrrction, only the lumped weights are 
considered as design variables. The results of optimization are snown in 
Table 14. The total weight Increases from 345 lb to 372 lb. The third 
collective mode of flapping and the second cyclic mode of inplane move 
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Table 14 Initial and final design for flapping, inplane and torsional mortis 
or Teetering rotor (Hltb pretnistl blade with boa bean dimension fixed 
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7. ARTICULATED ROTORS 

7.1 Definition 

Articulated rotor blade will be the subject of design blade in this 
section. Ref[25] The prlaary design variables are saae as before: the 
nail thickness of box beaa and luaped weights. Because the blade is 
articulated with a rigid hub. there is no distinction between collective 
or cyclic nodes for flapping and inplane. The blade is pretwisted. The 
boundary condition for flapping is a hinge at the root. There is root 
spring for torsional notions and an offset for inplane. Table 15 gives 
data for both the initial (and nlnlnun) blade stiffnesses and inertias 
as well as for the initial variables such as box beaa dinenslons, luaped 
weights. etc. 

7.2.1 Variable Box Diaenslon 

Box bean dinenslons. luaped weights and root spring are taken as 
design variables. Tables 16 shows the istial and final results of the 
optlnization procedure. The inplane frequency noves fron 4.84 to 4.69 
no/rev while the torsional node noves fron 4.25 to 4.45. The total 
weight drops slightly fron 96.55 to 95.48 lb. The root spring changes 
from 2. 41x10 ^ to 2.89x10*. 

7.3. Fixed Box Beaa Dinenslons 

Box beaa dinenslons will be considered to be fixed in this section. 
Only luaped weights and the root spring are taken as design variables. 
The final results are shown in Table 17. 



Rntatina Speed ; 293 RPM * 2 ' aa »9 of aoMent j 

Vounqs Modulus » 0.160 x 10 lb/ In inertia > 0.2000x10 lb-ln 

Axial Stiess < 30,000 psl density of box beaa i 0.000417 augs/ln 
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Table 17 Initial and final design for flapping, lnplane and torsional modes 
of Articulated rotor blade with box beam dimension fixed 
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Weighting Factor * * * * * * 1.00 

Desired Frequency A A * * 4.50 


8. RELATION BETWEEN VIBRAIICN AMD roPQTJFMCY 


8.1 FORMULATION 

In this section, we would like to show whether or not the forced 
response of the blade can he adequately controlled, as we have assumed, 
by our approach of 'frequency placement', that is, of restricting the 
natural frequencies of the blade to lie within narrow intervals located 
away fro* certain irteger multiple* of the rotor speed. Also we would 
examine whether or not aerodynamic damping substantially reducej the 
resonant peaks, in which case concern about avoiding resonances through 
proper selection of frequency windows would be unnecessary. Finally, the 
sensitivity of the optimal design to the choice of frequency vindow 
will be studied. 

This investigation is carried out through two, somewhat overlap- 
ping, problems. First, the forced response of an initial (i.e., non- 
optimized) design is eompar . to the response of a final design; cases 
with and without aerodynamic damping are considered. Next, the response 
of initial and final designs are evaluated as a single natural frequency 
is varied (the others being held fixed). In each ease, a forcing 
function containing harmonics of the rotor speed is applied. Again, 
cases with and without aerodynamic damping were considered. The general 
finding *rom these studies is that frequency placement is a viable means 
of reducing vibration, although it is by no means the only method and 
should be used in conjunction with others. 

8.2 RESPONSE VERSUS FREQUENCY 

The equations of motion for the finite-element representation of a 
rotor blade, subjected to an external excitation, may be written in 
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matrix fora as 

[MHX(t)> ♦ lcKX(t)} ♦ [XHX(t)> - {F(t)}, (27) 

where 

(Ml - wass matrix, 

<X(t)> ■ eo loan rector of nodal displacenant, 

[C] » dawping matrix, 

[X] - stiffness matrix, and 

{f(t)l ■ forcing function coltnn rector. 

The forcing function may in turn be expressed as 

<«t)> - {V q > e iwt , ( 28 ) 

where 

w -■ forcing frequency, and 
▼ Q ■ forcing amplitude. 

After some calculation, it can be shown that tbe amplitude of tbe 
response — written as <X>, independent of t — can be given as 

{X> - [X ♦ iwC - w 2 Hr l {V > ( 29 ) 

o 

In this section, the response of both the initial and final 

(optimal) designs to an external forcing function is studied as the 

frequency of tbe forcing function is raried. Blades both with and with- 
out aerodynamic damping are considered. To formulate these problems, 

consider tbe forced behavior of a rotor-blade. Only the flapping res- 
ponse la considered. The Inplane response is inferred from the results 


without damping, since there Is little aerodynamic damping in the in- 
plane direction. 

Fig. [13] shows a plot of the forcing amplitude V [Ref. 25] used 

o 

in the study. Given the forcing amplitude, we can calculate the response 

of each node of the finite-element representation of the blade as the 

value of the forcing frequency, w, is varied. The tip ( finite-element 

node farthest from the hub) response is of special interest. Before the 

results obtained from this study are presented, it is useful to examine 

the frequency placement results which are described in chapter 6 (see 

Table 10). The results for the frequencies (in units of cycles/rev) are, 
for flapping mode only. 


MODE 

INITIAL DESIGN 

FINAL DESIGN 

1st 

1.18 

1.18 

2nd 

3.22 

3.09 

3rd 

5.89 

5.67 


Blade dimensions are given in Table 9. 

The frequencies in the above table correspond to the symmetric 
modes of a teetering rotor. Thus, only even harmonics of the rotor speed 
have been considered as forcing frequencies. As a result, the optimized 
blade (Final design) finds the third mode moved away from the critical 
6.0/rev (from 5.89 to 5.67). Similarly, the movement of the second mode 
to 3.09/rev removes it from 2.0 and 4.0/rev. In the comparison study to 
follow, however, we will apply the entire spectrum of frequencies to 
th<s blade (not Just even harmonics). Thus, the 'Final Design' can no 
longer be considered optimum. A comparison of the two blades, however. 
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does indicate th« strong affect of resonance beoause each oase has a 
distinct resonanoe (6 and 3/rev). 

We shall now consider the results of the present study. Fig. 14 
shows the tip responses of both the Initial design and final (i.e. . 
optimized) design as functions of the forcing frequency. Aerodynes '■ 
damping has been negleoted (Alternatively, the results can be interpret- 
ed as giving the Inplane response). It can be seen that near 1.18 
eycles/rev, the responses of the two designs are very similar. However, 
the responses corresponding to the second and third modes differ signi- 
ficantly. For example, in the second mode, the peak of 3.22 cycles/rev 
(initial design) moves to 3.09 cycles/rev (final design)/ Similarly, the 
peak of the third mode moves from 5.89 oycles/rev to 5.(7 cycle/rev, 
which is especially important sinoe it is highly desirable to keep the 
frequency away from the integer frequency of 6 cycles /rev. We conclude 
from these results that the frequency placement approach does have a 
significant effect on the forced tip response when damping is not consi- 
dered. 

Next, the effect of aerodynamic damping is considered. The presence 
of the damping implies that results to be presented aorr*,jpond to flap- 
ping. Kathematioal details of the damping formulation are available 
[Ref. 181. The effect of aerodynamic damping on reducing the resonant 
peaks of the tip response of the initial blades is shown in Fig. 15. 
Fig. 16 shows the damped responses of both the initial and optimized 
blades so that the effect of frequency placement can be studied. It is 
interesting to observe here that when damping is included, no apparent 
advantage is gained by optimizing the blade, at least in terms of 
reducing the tip response, oxoept in the range of 3-4/rev, in which a , 
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Figure 15 Tip Response Versus Forcing Frequency for Initial 
Design Both With 2 nd Without Darrping 
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which a thirty-five percent reduction occurs. However, we oust also 
examine the effect of optimization when the response is measured by the 
average shear force existing in the blade. 

Consequently, the shearing force in the blade is considered next. 
As a measure of the average shear in the rotor, we consider the sum of 
the squares of the shear force (abbreviated SSS) , 


SSS 


V ♦ T 2 J * 


♦ Y 


10 


( 30 ) 


In this section, Y^ represents the shear force at node i in the (ten- 
element) finite-element model. "ate that the root shear is necessarily 
included as one of the terms on the right-hand side of the equation, so 
that a large value of root shear will cause SSS to also be large. 

Pig. 17 shows the variation of SSS with respect to the forcing 
frequency for the initial design with and without aerodynamic damping. 
Fig. 18 shows the same quantities for the final (optimized) design. Fig. 
19 cosipares the quantity SSS corresponding to initial and final designs 
when aerodynamics is considered. Inspection of these figures shows that, 
in contrast to behavior of the tip response, the shear response is 
significantly affected by changing blade frequency, even when 
aerodynamic damping is included. Thi 3/rev loads are increased by fifty 
percent due to the movement of v 2 from 3.22 to 3. 09/rev. Similarly, the 
6/rev loads are reduced by seventy percent due to the movement of w^ 

from 3.89 to 5.67/rev. Thus, even with damping, frequency placement is a 
powerful driver of loads. It follows that frequency placement can be 
justifiably considered an important part of blade optimization. 
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Figure 17 Sub of Squares of Shears Versus Forcing Frequency 

for initial Design Botn With and Witnuut Damping 
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8.3 RESPONSE VERSUS PLACEMENT 

In the study just described, the response of the blsde to changes 
in the forcing frequency was considered. Nov we consider s different 
approach. In effect, ve examine how the blade responds to a forcing 
function 'during the optimization procedure' — in the sense that during 
optimization, the optimization algorithm varies the natural frequency of 
the blade (to force it to satisfy the frequency constraints). In obtain- 
ing the results to be presented next, ve simulated the optimization 
procedure by varying the natural frequency. Thus ve can observe vhat 
happens to the forced response during frequency placement. 

The formulation of the approach is as follovs. Through appropriate 
transformations (described in Appendix 11.3), the system mass matrix can 
be written as 

[mj - lurtur 1 (3D 

and the system stiffness matrix as 

IK] - iMjfUl diagl (v. 2 ) )tU] T [M] ( 32 ) 
were v^ are the natural frequencies of the system, [U] is a matrix whose 
columns are eigenvectors, and the notation 'diag' indicates a diagonal 
matrix (all off-diagonal terms vanish). From examination of these 
expressions, it can be seen that the stiffness and mass matrices can be 
considered functions of the natural frequencies. Thus it becomes 
possible to fix all frequencies but one, and then study the response of 
the system as that one frequency is varied with mode shapes also held 
fixed. In particular, the response to the following forcing function 


14 * 


will be studied: 


<F(t)> - {V^ iwt ♦ {7 2 }, i2jrt 


( 33 ) 


♦ ... ♦ {Vje^ 

vhere 


v * the rotor speed, and 
{▼,} - (I/nHV }. 

D O 


( 34 ) 


where {V q } was defined previously in Fig. 13. Since the arguments of the 
exponentials are integer multiples of w, resonance will occur at har- 
monics of the rotor speed. The particular forcing function given above 
is known from empirical observation to provide an approximate, but 
physically realistic representation of the radial and harmonic variation 
of the amplitude of the load on a real blade. As mentioned in this 
section, the blade response will be defined by the tip displacement and 
the sum of the squares of the shears except that, here, the n * 1 term 
has been omitted from the expressions for calculating tip displacement 
and shears because this term represents a tip-path plane tilt that is 
controlled by the pilot fc:- trimming purposes. It is not part of the 
true vibratory loads we are considering. 

Results for the problem just formulated are shown in Fig. 20, where 
the sum of the squares of the shears is plotted as a function of w^, the 
second natural frequency, with the other natural frequencies being 
fixed. This figure corresponds to the Initial blade design (blade dimen- 
sions are given in Table 9. Fig. 21 shows the same quantity for the 
case where the third natural frequency is varied. It Is interesting to 
note that the response curve '.or the damped case in Fig. 21 lacks 
resonant peaks — apparently '.he damped response is so completely 
dominated by the resonance of the second natural frequency, which is 
fixed near 3/rev, that the (damped) resonant peaks for the third fre- 


quency are negligible by comparison. It should be noted that the 
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Figure 20 Sue of Squares of Snears Versus Second Natural Frequency 
for Initial Design Both With and without Dadoing 
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response of the original design Is respresented by only one point on 
Figure 20 or Figure 21. This value can be found by taking w a ■ 3.22 or 
w ? - 5.89 on those figures, sinoe these are the initial design values 
(Table 10). 

For the final (optimal) design, the analogous quantities are plot- 
ted in Figs. 22 and 23. Again, no resonant peaks are present in the 
damped response when the third natural frequency is varied. Comparison 
of magnitudes of ordinates in Figs. 20 and 22 (no damping) shows that 
the overall shear measure is reduced in the final design in the regions 
away from resonance. Also, the choice of scale or. the vertical axis in 
Fig. 22 highlights the effect of frequency placement. Note that by 
inspection of Figs. 20-23. a designer may select the design frequency 
which minimizes the average shear as measured by the SSS. 

One of the most interesting results of Fig. 22 is information about 
the width of valleys and peaks, since this gives design information. 
First, let us examine the non-damping curve (inplane response). Here, 
the minimum points are nearly at the centers of the regions (2.55/rev) 
and (3.5 f /rev). The frequency windows to maintain no more than thirty 
percent increase in loads are 2.40 - 2.70/rev and 3.40 - 3.70/rev (plus 
or minus 0.15/rev) — a fairly narrow window. For the damped curves 
(flapping response), mlmima are also near the one-half points, but the 
window for thirty-percent increases are much wiuer: 2.20 to 2.90/rev and 
3.20 to 3.80/rev (plus or minus 0.30/rev). Stated another way, inplane 
frequencies should be no closer than a 0.4/rev from integers, but flap- 
ping frequencies may be as close as 0.2 from an integer. It should be 
emphasized that these observations apply to this particular example and 
a, ay not be generalized for other frequency constraints. 
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Figure 23 Sum cf Squares of Shears Versus Third Natural 
for Finji Design Born Wien and Wirnour D.-mpm 
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Another conclusion to be drawn from the above results is that the 
undamped response curve has very flat-bottomed 'valley' when one of the 
fixed frequencies is near an integer value (of. Figs. 22 and 23). 

8.4 RESPONSE DDE TO EVEN-INTEGER HARMONICS 

In this section, we still study the response of the blade due to 
the ohange of forcing frequency. However, the forcing frequency will be 
a little different from the one in section 8.3.1* We will consider only 
the even Integer forcing frequencies for which the collective me are 
optimized. Thus, the forcing funotion may be written as 

(P(t)l - V 2 . iJ,rt * » 4 . i4,rt * T s . 1Swt * T,. 1 *"* 
where 

w - rotor speed 

V n - 1/n • V 
n o 

Since the arguments of the exponentials are even integer multiples of w 
only, the resonance will expectedly occur at an even integer of 
harmonics of the rotor speed. The response of shear stress will be 
studied in this section. 

Results for this case are shown in Fig 24, where the sum of the 
squares of the shear harmonics is plotted as a function of w 2 (w 2 and Wj 

are fixed; Wj ■ 1.18 and w^ - 5.89). The figure corresponds lo the 
initial blade design. Resonance peaks oocur at 2 and 4/rev. Thus, for 
the shear response, a second frequencies around 3.0 would be the 
suitable choioe for w^ %s design frequncy. Fig 25 is a similar 
comparison for a variable value of Wj ( w^ - 1.18 and w 2 - 3.22 are 
fixed). Here, the resonanoe is at 6.0/rev, and Wj ■ 5 or Wj ■ 7 would be 
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Figure 24 


Sum of Squares of Shears Versus Second Natural Frequen> 
for Initial Design Both With and Without Damping 
(where Forcing Function is Even Integer Multiple) 
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fWITl^L DESI3N fUITH AND WITHOUT A£PO OVN*M I C S ) 



Figure 25 Sum of Squares of Shears Versus Third Natural Frequency 
for Initial Des' Both with and without Damping 
(where Forcing . -.action is Even Integer Multiple) 
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Results are shown in Fig 26 and 27 for the final design blade. 
Again, it is noted that the resonance peaks appear at even integers. A 
comparison of Fig 24 and 25 with Fig 26 and 27 shows the relatively 
lower vibrations of the final design. In either case, however, we see 
the sensitivity of vibrations to frequency, even with aerodynamics. 
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9. SUMMARY AND CONCLUSION 

The optimization technique works very successfully on the design 
of rotor blades even when there are as many as 55 constraints. The most 
efficient optimization orocedure involves 2 steps. In the first step, 
the objective function is based on frequency placement with appropriate 
structural constraints. In the second step, the objective function is 
weight with frequency windows as constraints. 

A? far as the optimization of helicopter blades is concerned, the 
appropriate constraints include autorotatlonal inertia, axial stress, 
geometric limitations of the cross-section, and the placement of mass 
center 'I'orwaid of the quarter chord. 

Proper choice of input data can ensure the optimization runs 
smoothly and converges faster. Although we have up to 15 full con- 
straints and 40 side constraints at one time, the program works very 
well. The reason may be due to the fact that the input data are 
practical enough to meet (or to be close to) most of the constraints 
even before the optimization starts. However, if we start the optimiza- 
tion with random input data, the results may not be as good as expected. 

The forced response of the blade can be adequately controlled, as 
expected, by the approach of 'frequency placement'. 

The optimization techniques results in realistic designs by place- 
ment mass at antinodes or nodes, by adding stiffness at antinodes or 
nodes, and by placing mass near tne tip to achieve autorotational 
inertia at minimum weight. 
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APPENDIX 11.1 

Nonane lature 

size of tip mass 

area of eroas-tection of box bean 

width of box beam 

blade chord 

damping matrix 

width of lumped mass 

wall thickness of box beam 

weighting function 

flapping stiffness 

inpiane stiffness 

safety factor 

final design 

forcing function 

gravity 

torsional stiffness 
height of box bean; 
initial design 

area moment of inertia of box beam 

constant area moment of inertia of blade 

mass moment of inertia of box beam 

constant mas moment of inertia of blade 

stiffness matrix 

area moment of inertia 

mass moment of inertia 










1 

[Ml 


pi ,p2 ,p3 


r 

R 


3 

sss 

t 



0 





v 


W 1 ,w n 
wf 

x 

e 

<p 

/* 

e 

V 


length of an element 
mas 8 matrix 

blade first three frequencies (no/rev) 

distance from rotation axis 

length of the blade 

box beam thickness 

sum of squares of shear force 

box beam thickness 

tension force 

kinetic energy 

total energy 

displacement of degree of freedom of element 

unit matrix 

forcing amplitude 

displacement in yx plane' 

displacement in zx plane 

blade frequencies 

weighting factor 

coordinate axis and length parameter 

blade twist 

pretwist angle 

mass density of the blade 

small parameter 

nt 


fl rotating speed 

maximum stre'"' 
C mass of skin 
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density of lumped mass, box beam, honeycomb 
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APPEUDIX 11 .1 

Calculation of Torsional Stiffness, GJ 
The blade cross section is shown in Fig. 1 and is idealized into a 
two cell torsion box in Fig. 28. Although the bending and torsional 
Inertias inolude the contribution of all masses (including the filler 

elements), the torsional stiffness is based only on the structural box 

and thln-slcin elements. The torsional stiffness, GJ, is based on classi- 
cal thin-walled closed section theory. (Ref .26) The effects of warping re 
and distortion of the cross section are neglected. The equations are 
briefly summarized as follows: 

Considered cell i, having an enclosed area of A^, and thickness t. 
The length along the circumference is measured by c, and the shear 
stress is denoted by s.. The torsion oonstant is 

J - [2A / f (s dc)]T ( 37 ) 

1 * I 

where J- t do is the line integral along the entire closed box, and T 

represents the total torque applied to the entire system. We also have 

2TA /J - f (s do) - f (q dc/t) 
l J) J \ 


( 38 } 


where c is the constant shear flow in the i-th box. Denoting the single 
adajcent box as k, there will be a common wall, ik, between the two 
boxes. The shear flow in the coenon wall will be the difference of shear 
flows q^ and q fc . Thus, for the i-th box, 

q. f (dc/t) - q. f (dc/t) - 2TA./J (39) 

1 i K ik 1 

and for the k-th box 

-’i * ’1 J k (ic/t) ■ “n''' (40, 

Defining q. ■ Q.(2T/J) and n. “ f (dc/t) , we can reform 
J J 

Eqs.3 to read 

n ii Q i “ n ik\ “ A i (41a) 

“°ik Q i + ^k^ * \ (4lb) 

Furthermore , 

T - 2(q.A. + q k A fc ) - (4T/J)(Q.A. ♦ Q^) (42) 

Thus, the torsional constant, J, is given by 


4(Q.A. 
i l 


W 


(43) 


Specifically, for the two cell box shown in Fig. 28, cal ber 1 the 
box beam and box 2 the equivalent trailing edge. Denote the 
circumference ot box 2 by C. The areas, Aj and A 2 


are 


given. 



Then, from the other dimensions 


fljj ■ n/dj + 2(b/tj) * h/d 2 

(44a) 

n 12 “ 

(44b) 

n 22 - C/p ♦ h/d 2 

(44c) 

Then, solving Eqs. 4e tnd 4b gives 


Qj - l Aj Q 22 + *2 n 12^ D 

(45a) 

Q 2 ■ (A 2 n^ + A^nj 2 )/D 

(45b) 

D ■ B U a 22 * ( "12 )2 

(45c) 


That, til qatntiti.es in Eq.6 tre known tnd the tortiontl const tot, J, 


ctn be evtlutted 


Cell I: Constant- thicknesses for sides of box beam 
represent weighted effects of variable 
thickness elements. 

Cell 2 : This cell represents configuration cf 
skin and trailing edge. 


Idealized tvo-cell model for calculation of torsional 


stiffness . 
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APPENDIX 11.3 

Derivation of Mas* and Stiffness Matrices 
aa Functions of Natural Frequencies 







Def ine 

[K*I - [M]“ 1/2 [K][M]' 1/2 , ( 46 ) 

and construct a square matrix [D ] by using the eigenvectors of [K ] as 
columns. If the eigenvectors are normalized to the identity matrix, that 
is, if 

[U*] T tO*] - [I], (4 7) 


i„ then follow* that 


tU*] T [K*][U*l - diag i (w. 2 ) 1, 


where v. 

i 


Next, 


are the eigenvalues of [X ] . 
let 


[or] - [m]" 1/2 [d*]. 


from which it follows that 

lu! T - [aVlM]' 172 , 

Eul” 1 - [u*] T [M] 1/2 , 

fn] T (Kl[0] - diag l (w. 2 ) 1 , 
end 

[U] T [M] fU] - [Ij. 


( 48 ) 


( 49 ) 


( 50 ) 


( 51 ) 

( 52 ) 


( 53 ) 
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Finally., then, the stiffness and mess matrices can be written as 

2 

functions of the eigenvalues, w^ : 

ini - £0 ; _t [oj~ 1 

- !54 ’ 

and 

[K] - [ur T diag[ (w. 2 ) Hu]’ 1 

- [M] 1/2 [0*]diag[ (w. 2 ) ][0*] T [M] 1/2 

- [M] [UJdi.tgt (w. 2 ) ][U] T [M]. ( 55 ) 

•J 

Note that the eigenvectors, [U] , and eigenvalues, , appearing on 
the right-hand side were originally calculated from the stiffness and 
mass matrices, [K] and [Ml . If we consider only relatively small changes 
in the frequencies, w^, then the eigenvectors should relatively 
unchanged. Thus the last two equations for iM] and [K] with [U] held 
fixed can be considered as expressing the mass and stiffness matrices as 
explicit functions of the natural frequencies. 
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